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Instead of the fermion-fermion scattering which identifies the BCS-BEC crossover in cold atom 
systems, boson-boson scattering is measurable and characterizes the BCS-BEC crossover in QCD 
superfluid. We study n-n scattering in a pion superfluid described by the Nambu-Jona-Lasinio 
model. We found that the amplitude of the scattering length drops down monotonously with 
decreasing isospin density and finally vanishes at the boundary of the phase transition. This indicates 
clearly a BCS-BEC crossover in the pion superfluid. 
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There are two kinds of condensed states in usual 
fermion gas, the Bardeen-Cooper-Shrieffer condensa- 
tion (BCS) of fcrmions where the pair size is large and 
the pairs overlap each other, and the Bose-Einstein 
condensation (BEC) of molecules where the pair size 
is small and the pairs arc distinguishable. The BCS 
wave function can be generalized to arbitrary attraction 
which leads to a smooth crossover from BCS to BEC [l| . 
In cold atom systems, the experimental observable to 
identify the BCS-BEC crossover is the s-wave scatter- 
ing length between fermions, which is positive in BEC 
state and negative in BCS state 0. 

Recently the study on quantum chromodynamics 
(QCD) phase structure is extended to finite isospin 
density. For a QCD system at finite temperature and 
baryon and isospin density, the phase transitions in- 
clude not only color deconfinement Q, chiral symme- 
try restoration 0] and color superconductor Q, but 
also pion superfluid @. The increasing isospin density 
induces a phase transition from normal nuclear matter 
to pion superfluid, due to the spontaneous isospin sym- 
metry breaking. The critical isospin chemical poten- 
tial at zero temperature and baryon density is exactly 
the pion mass [fj. By analogy with the usual super- 
fluid, the BCS-BEC crossover in pion superfluid can be 
theoretically described @, H| by the quark chemical po- 
tential which is positive in BCS and negative in BEC, 
the size of the Cooper pair which is large in BCS and 
small in BEC, and the scaled pion condensate which 
is small in BCS and large in BEC. However, unlike 
the fermion-fermion scattering in cold atom systems, 
quarks are unobservable degrees of freedom, and thus 
the quark-quark scattering can not be measured or used 
to experimentally identify the BCS-BEC crossover. 

In pion superfluid, the pairs themselves, namely the 
pion mesons, are observable objects. One can measure 
the 7r — it scattering to probe the properties of the pion 
condensate and in turn the BCS-BEC crossover. Since 
pions are Goldstone modes corresponding to the chiral 
symmetry spontaneous breaking, the it — 7r scattering 
provides a direct way to link chiral theories and experi- 
mental data and has been widely studied in many chiral 
models 0, [13] • Note that pions are also the Goldstone 
modes of the isospin symmetry spontaneous breaking, 
the 7r — 7T scattering should be a sensitive signature of 
the pion superfluid phase transition. 

While the perturbative QCD can well describe the 
properties of the new phases at extremely high temper- 
ature and density, the study on the phase transitions at 



moderate temperature and density depends on lattice 
QCD calculations and effective models with QCD 
symmetries. One of the widely used effective models 
is the Nambu-Jona-Lasinio (NJL) model [Hj], which is 
originally inspired by the BCS theory and its version 
at quark level gives simple and direct description 
of the dynamic mechanisms of spontaneous chiral sym- 
metry breaking, color symmetry breaking and isospin 
symmetry breaking. The s-wave tt — tt scattering length 
calculated llCj in the model is consistent with the Wein- 
berg limit [14| and the experimental data fl5j in vac- 
uum. In this Letter, we extend the calculation to finite 
isospin chemical potential and focus on its relation to 
the BCS-BEC crossover in the pion superfluid. 

The Lagrangian density of the two flavor NJL model 
at quark level is defined as [l3| 

C = V> {vy^d^ - m + jofi)ip + G (ipip) 2 + (tpi^nip) 

(1) 

with scalar and pseudoscalar interactions correspond- 
ing to a and 7r excitations, where mo is the current 
quark mass, G is the four-quark coupling constant 
with dimension GeV -2 , Tj (i = 1,2,3) are the Pauli 
matrices in flavor space, and fi = diag (fi u , fid) = 
diag (/is/3 + fii/2, /is/3 — fii/2) is the quark chem- 
ical potential matrix with fi u and fid being the u- 
and (i-quark chemical potentials and \ib and fii the 
baryon and isospin chemical potentials. At fii = 
0, the Lagrangian density has the symmetry of 
J7s(l) SUi(2) SUa(2), corresponding to baryon, 
isospin and chiral symmetry. At fij ^ 0, the symme- 
tries SUi(2) and SUa{2) are firstly explicitly broken 
down to Ui(l) and {/^(l), and then the nonzero pion 
condensate leads to a spontaneous breaking of Ut(1), 
with pions as the corresponding Goldstone modes. At 
fiB =0, the Fermi surface of u(d) and anti-d(u) quarks 
coincide and hence the condensate of u and anti-ci is 
favored at /.*/ > and the condensate of d and anti-it 
quarks is favored at fii < 0. Finite fiB provides a mis- 
match between the two Fermi surfaces and will reduce 
the pion condensation. 

Introducing the chiral and pion condensates a = 
(il^ip) and 7r = (ifti'ysTi'ip) and taking them to be real, 
the quark propagator S in mean field approximation 
can be expressed as a matrix in the flavor space 
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with the dynamical quark mass m q = mo — 2Go gen- 
erated by the chiral symmetry breaking. By diagonal- 
izing the propagator, the thermodynamic potential can 
be simply expressed as a condensation part plus a sum- 
mation part of four quasiparticle contributions @ . The 
gap equations to determine the condensates a (or quark 
mass m q ) and it can be obtained by the minimum of 
the thermodynamic potential. 

In the NJL model, the meson modes are regarded as 
quantum fluctuations above the mean field. The two 
quark scattering via meson exchange can be effectively 
expressed at quark level in terms of quark bubble sum- 
mation in RPA fijjj. The quark bubbles are defined 
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with indexes m,n — a, tt+ , 7r_ , ttq , where the trace 
Tr = TrcTrpTrjj is taken in color, flavor and Dirac 
spaces, the four momentum integration is defined as 
fd i p/(2n) i = iTYsj fd 3 p/{2n) 3 with fermion fre- 
quency po = iujj = i(2j + 1)ttT (j = 0, ±1, ±2, • • •) at 
finite temperature T, and the meson vertices are from 
the Lagrangian density (JTJ , 



1 m = a 
. ij 5 r+ m = 7T+ 

x *75 r - m = 7r - 
m = 7T , 



r* = 



1 m = a 
i"fsT— m = 7r+ 

i757"3 m = 7T 



(4) 

Since the quark propagator S contains off-diagonal el- 
ements, we must consider all possible channels in the 
bubble summation in RPA. Using matrix notation for 
the meson polarization function n(fc) in the 4x4 meson 
space, the meson propagator can be expressed as 
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Since the isospin symmetry is spontaneously bro- 
ken in the pion superfluid, the original meson modes 
a, 7r + , 7r_ , 7To with definite isospin quantum number 
are no longer the eigen modes of the Hamiltonian 
of the system, the new eigen modes ef, 7F+ , tF - , ttq 
are linear combinations of the old ones, their masses 
Mi(i — £j, 7T_j_ , 7T_ , 7T0) are determined by poles of 
the meson propagator at kg — Mi and k = 0, 
det [1 — 2GH(Mi, 0)] = 0, and their coupling constants 
gi q q are defined as the residues of the propagator at the 
poles. 

The condition for a meson to decay into a q and a q 
is that its mass lies above the q — q threshold. From the 
pole equation, the heaviest mode in the pion superfluid 
is a and its mass is beyond the threshold value. As a 
result, there exists no er meson in the pion superfluid, 
and the coupling constant g-^qq drops down to zero at 
the critical point /ij and keeps zero at uj > (l6j . 

We now study ir — ir scattering in the pion super- 
fluid. To the lowest order in 1/N C expansion, where N c 
is the number of colors, the invariant amplitude T of 
the it — ir scattering is calculated from the box diagrams 
in s, t and u channels, as shown in FigUfor the s chan- 
nel. Different from the calculation in normal state fltl ] , 
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FIG. 1: The lowest order diagram for n — n scattering in the 
pion superfluid. The solid and dashed lines are respectively 
quarks and pions, and the dots denote meson-quark vertices. 



the (T-exchange diagrams vanish here due to the dis- 
appearance of the a meson in the pion superfluid. In 
the limit of the scattering at threshold yfs = 2M W and 
t = u = 0, the amplitude T approaches to the scatter- 
ing length a (jn^ J ) in unit of the inverse of dynamical 
pion mass through the relation 
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Note that the threshold condition can be fulfilled by 
a simple choice of the pion momenta, k a = kb = 
k c = kd = k and k 2 = At£ = s/4, which facili- 
tates a straightforward computation of the diagrams. 
Here the meson dynamical mass is obtained from 
the quark and meson pole masses and the quark dy- 
namical mass. For instance, we have niw + = My, + 

2 (m, - v /K-w/2) 2 +4GV). 

In the calculation in normal matter, people are inter- 
ested in the 7r-7r scattering length with definite isospin, 
a i=o,2, which can be measured in experiments due 
to isospin symmetry. However, finite isospin density 
breaks down the isospin symmetry and makes the scat- 
tering length ao,2 not well defined. In fact, the new 
meson modes in the pion superfluid do not carry def- 
inite isospin quantum numbers. Unlike the chiral dy- 
namics in normal matter, where the three degenerated 
pions arc all the Goldstonc modes corresponding to 
the chiral symmetry spontaneous breaking, the pion 
mass splitting at finite /zj results in only one Gold- 
stone mode 7f+ with M w+ = in the pion superfluid. 
In the following, we focus on the scattering process of 
7F+ + 7T+ — > 7r+ + 7T_|_ and calculate the scattering 
length a+. 

The scattering amplitude for any channel of the box 
diagrams can be expressed as 
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with the quark propagators Si = S3 = S(p), S2 = 
S{p+k), and S4 = S(p—k) for the s and t channels and 
Si = S3 = S(p + k) and S2 = S4 = S(p) for the u chan- 
nel. In the limit of threshold, k = (M^ + , 0) = (0, 0) in 
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the pion superfluid leads to S\ = S2 = S3 = S± = S(p) 
for all the three channels. Making the fermion fre- 
quency summation over the internal quark lines, the 
total amplitude at the threshold is simplified as 
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where E^. = E± =F //b/3 are the energies of the four 



quasiparticles with E± = J (E ± fii/2) + 4G 2 tt 2 ,E = 

p 2 +mg, and f(x) is the Fermi-Dirac distribution 

function f(x) = (e x l T + 1) 1 with the derivative 
f(x)=df/dx. 

Since the NJL model is non-renormalizable, we can 
employ a hard three momentum cutoff A to regular- 
ize the gap equations for quarks and pole equations 
for mesons. In the following numerical calculations, we 
take the current quark mass mo — 5 MeV, the cou- 
pling constant G = 4.93 GeV -2 and the cutoff A = 653 
MeV[l7|. This group of parameters correspond to the 
pion mass m v = 134 MeV, the pion decay constant 
/ w = 93 MeV and the effective quark mass M q = 310 
MeV in the vacuum. 
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FIG. 2: The scattering length \a+\ for the process 7F+ + 
7f + — > 7f+ + n+ as a function of isospin chemical potential 
\ii in the pion superfluid phase. 

In Fig 121 we plot the scattering length \a+\ as a func- 
tion of isospin chemical potential [i\ in the pion super- 
fluid phase at two temperatures T = and T = 100 
MeV, keeping baryon chemical potential fis — 0. From 
(O , the scattering length is governed by the meson cou- 
pling constant, a+ ~ 7+ ~ 9^ + qq- ^ n * ne pi° n super- 
fluid, the meson mode W + is always a bound state, its 
coupling to quarks goes up with increasing fii [lj| , and 
therefore the amplitude of the scattering length \a + \ 
drops down with decreasing fii and reaches zero at the 
critical value /xj due to ffw + qq — at this point, where 
M J = = 134 MeV at T = and fij = 142 MeV at 
T = 100 MeV. 



The above ^/-dependence of the meson-meson scat- 
tering length in the pion superfluid can be understood 
well from the point of view of BCS-BEC crossover. We 
recall that the BCS and BEC states are defined in the 
sense of the degree of overlapping among the pair wave 
functions. The large pairs in BCS state overlap each 
other , and the small pairs in BEC state are individual 
objects. Therefore, the cross section between two pairs 
should be large in the BCS state and approach zero in 
the limit of BEC. From our calculation shown in FigJ5J 
the 7T — 7T scattering length is a characteristic quantity 
for the BCS-BEC crossover in pion superfluid. The 
overlapped quark-antiquark pairs in the BCS state at 
high isospin density have large scattering length, while 
in the BEC state at low isospin density with separa- 
ble pairs, the scattering length becomes small. This 
provides an experimental observable for the BCS-BEC 
crossover at QCD level, analogous to the fermion scat- 
tering length in cold atom systems. 
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FIG. 3: The scattering length |a+| for the process 7f+ + 
7f+ — » ¥+ + 7f+ as a function of temperature T in the pion 
superfluid phase. 

With increasing temperature, the pairs will gradu- 
ally melt and the coupling constant g^qq drops down in 
the hot medium, leading to a smaller scattering length 
at T — 100 MeV in comparison with the case at T = 0, 
as shown in FigJ5J To see the continuous tempera- 
ture effect on the scattering length in the BCS and 
BEC states, we plot in Fig|3] \a + \ as a function of T 
at fii = 160 and /i/ = 400 MeV, still keeping (j,b = 0. 
While the temperature dependence is similar in both 
cases, the involved physics is different. In the BCS 
state at [ij = 400 MeV, \a+\ is large and drops down 
with increasing temperature and finally vanishes at the 
critical temperature T c = 188 MeV. Above T c the sys- 
tem becomes a fermion gas with weak coupling and 
without any pair. In the BEC state at fj,j = 160 MeV, 
the scattering length is much smaller. At a lower crit- 
ical temperature T c = 136 MeV, the condensate melts 
but the still strong coupling between quarks makes the 
system be a gas of free pairs with \a+ \ =0. 

There are some other quantities to describe the BCS- 
BEC crossover in superfluid or superconductor Q, 
which arc difficult to be experimentally measured but 
can be used to confirm the BCS-BEC crossover picture 
characterized by the pion scattering length. We calcu- 
late the scaled binding energy e//x/ as a function of /ii 
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FIG. 4: The scaled 7r+ binding energy e/fj,i as a function of 
isospin chemical potential /i/ in the pion superfluid. 

in pion superfluid at T = and T = 100 MeV, shown 
in FigHJ The binding energy of 7f+ is defined as the 
the mass difference between 7f + and the two quarks, 
e = M w+ — M u — M-g. With decreasing isospin chemi- 
cal potential, the binding energy increases, indicating a 
BCS-BEC crossover in the pion superfluid. When the 
medium becomes hot, the condensate melts and the 
pairs are gradually dissociated. These results support 
our conclusion of BCS-BEC crossover identified by the 
measurable meson scattering length. 

In summary, we proposed the meson-meson scat- 
tering as an experimentally measurable probe of the 
BCS-BEC crossover at quark level. Different from the 
fcrmion-fcrmion scattering which is often used to ex- 



perimentally identify the BCS-BEC crossover in cold 
atom systems, quark scattering can not be measured 
and its function to characterize the BCS-BEC crossover 
in QCD superfluid is replaced by the meson scatter- 
ing. In the BCS quark superfluid, the large and over- 
lapped pairs lead to large pair-pair cross section, but 
the small and individual pairs in the BEC superfluid in- 
teract weakly. In the frame of a two flavor NJL model 
at finite temperature and isospin density, we calculated 
the 7r — 7r scattering length in the pion superfluid. It 
is large at high isospin chemical potential and drops 
down monotonously with decreasing isospin chemical 
potential and finally approaches zero at the border of 
the pion superfluid, indicating a BCS-BEC crossover. 

The meson scattering length \a+\ shown in Figs|2]and 
[3] are obtained in a particular model, the NJL model, 
which has proven rather reliable in the study on chi- 
ral, color and isospin condensates at low temperature. 
Since there is no confinement in the model, one may 
ask the question to what degree the conclusions can be 
trusted. From the general picture for BCS and BEC 
states, the feature that the meson scattering length ap- 
proaches to zero in the process of BCS-BEC crossover 
can be geometrically understood in terms of the degree 
of overlapping between the two pairs. Therefore, the 
qualitative conclusion of taking meson scattering as a 
probe of BCS-BEC crossover at quark level may survive 
any model dependence. 
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